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We analyse the origin of the magnetic field decay in normal radio pulsars found by us in a recent study. This decay has a 
typical time scale ~ 4 x 10® yrs, and operates in the range ~ 10® - fewx 10® yrs. We demonstrate that this held evolution 
may be either due to the Ohmic decay related to the scattering from phonons, or due to the Hall cascade which reaches 
the Hall attractor. According to our analysis the hrst possibility seems to be more reliable. So, we attribute the discovered 
held decay mainly to the Ohmic decay on phonons which is saturated at the age fewx 10® yrs, when a NS cools down to 
the critical temperature below which the phonon scattering does not contribute much to the resistivity of the crust. Some 
role of the Hall effect and attractor is not excluded, and will be analysed in our further studies. 
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1 Introduction 

Observational appearances of most (if not all) types of 
neutron stars (NSs) significantly depend on their magnetic 
fields. The magnetic field properties include not only its 
strengths, but also its topology. Evolution of both was 
a subject of numerous studies (see an early review in 
|Chanmugam 1992| and more recent in |Geppert 2009[ ). 

Theoretical models typically predict that the field might 
decay due to several reasons. However, analysis of obser¬ 
vational data mostly yields controversial results. Especially, 
if we restrict our consideration to the normal radio pulsars 
only. On one hand, decay of magnetar fields seems to be in¬ 
evitable for explanation of observational manifestations of 
these sources (see a review in Mereghetti, Pons & Melatos 
2015). Another strong evidence in favour of field decay 
comes from the fact that fields of millisecond (recycled) 
pulsars are much weaker than those of normal ones, which 
can be explained by significant decay due to accretion in 
low-mass X-ray binaries (LMXBs). On other hand, there is 
no definitive answer if standard magnetic fields (^ 10^^ G) 
significantly decay on the time scale few million years - the 
lifetime scale of a normal radio pulsar. 

Recently we proposed a modified data analysis to 
probe the magnetic field decay in normal radio pulsars 
( Igoshev & Popov 2014| l. Application of the new method 
resulted in discovery of magnetic field decay by a factor ^ 2 
in the time interval ^ 10^ - 3 x 10® yrs. It was speculated 
that on longer time scales this decay might stall in order to 
be in correspondence with previous studies which had not 
found any significant field decay on longer time scales. In 
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this paper we present preliminary results of our study on 
the reason for such non-uniform field evolution. The full 
consideration of the problem will be presented elsewhere 
(Igoshev & Popov, in prep.|l 


2 Field decay in radio pulsars 


In case of normal radio pulsars their magnetic field can 
not be directly measured. However, the basic pulsar prop¬ 
erties such as period and period derivative should already 
contain some information about magnetic field evolution. 
Still, there is no straightforward way to extract the field evo¬ 
lution because both the kinematic magnetic field estimate 


PB = \/Pp, where /3 is a (probably time dependent) coef¬ 
ficient determined by NS properties, and the spin-down age 
T = P/ (2P) are related to each other. It means that there is 
no good independent age estimate applicable to significant 
number of pulsars of different ages. Thus, it is necessary to 
use specific (statistical) methods to uncover magnetic field 
evolution from the data based on model-independent ob¬ 
servable parameters, such as P and P. 

In the previous study ( Igoshev & Popov 20I4| | we pre¬ 
sented a new path towards breaking the dependency be¬ 
tween the spin-down age and the field estimate. A new 
method to study field evolution in normal radio pulsars was 
proposed. The method is extensively tested with population 
synthesis technique using synthetic populations with differ¬ 
ent magnetic field evolution. 

Our method is similar to the well-known pulsar current 
approach (Vivekanand & Narayan 19811, but we study the 
“flow” of pulsars not along the direction of growing spin 
periods, but along the spin-down age. We propose to count 
a number of pulsars Ni in spin-down age intervals Ar^. If 
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Fig. 1 Magnetic field decay in normal radio pulsars from 
the Parkes Multibeam and Swinburne surveys (PMSS) ac¬ 
cording to Igoshev & Popov (2014). The obtained result 
cannot be directly extrapolated to smaller or larger ages. 


the birth rate is constant and selection effects are similar for 
pulsars with different spin-down ages, then the ratio of pul¬ 
sars number at different intervals of spin-down age refiects 
the ratio of true ages of pulsars in these intervals. 

With currently available data our method is applicable 
to sources with ages ^ (1-5) x 10® yrs. We obtained that 
the field decays in this age interval by a factor of ^ 2. 
This can be fitted by an exponential decay with time scale 
r « 4 X 10® yrs. We show the main result by Popov & 
Igoshev(2014) in Fig[T| There (and in other plots) The mag¬ 
netic field is represented as B{t) = Bq x f{t), where f{t) 
is the decay function which describes field evolution. 

Obviously, the result presented in Fig. [T] is not in corre¬ 
spondence with most of the previous studies (see, for exam¬ 


ple, Faucher-Giguere & Kaspi 2006 and references therein) 


if it is assumed that the field continues to diminish with the 
same rate at larger ages. On other hand, moderate field de¬ 
cay in normal pulsars is not in contradiction with observa¬ 
tions jGullon et al. 2014l l. So, the stage of relatively rapid 
decay might be terminated. The method is statistical in its 
nature, and therefore it cannot shed light on the physical 
mechanism which governs the magnetic field decay at this 
interval of true ages. In the new study we address this ques¬ 
tion. 

3 Hall cascade vs. Ohmic dissipation on 
phonons 


Basing on current understanding of NS crust properties, 
there are two effects which can cause the magnetic field de¬ 
cay with a time scale ~ 10® years: the Hall effect and Ohmic 


decay (Geppert2009 Gumming, Arras & Zweibel 2004). 


The former effect is not dissipative, it actually works simi¬ 
lar to a turbulent cascade and redistributes the magnetic field 
energy from large scales (for example, dipole component) to 
smaller scales (in particular, if initially there existed I multi¬ 
pole, then the attractor would consist of I and I -f 2 poloidal 


Fig. 2 Magnetic field decay due to scattering from impu¬ 
rities of crystalline lattice. The impurity factor is Q = 0.05. 


components, see |Gourgouliatos & Gumming 2014b| for de¬ 
tails). Intensity of the Hall effect depends on the value of 
initial magnetic field. 

Magnetic field decay in young NSs can be non¬ 
monotonic. There are several reasons for that. The first 
possibility is related to a specific property of the Hall 
cascade. The rate of field decay obtained in the arti¬ 
cle by [Igoshev & Popov 2014[ is close to the expected 
time scale of the Hall drift for normal magnetic fields 

was demonstrated by 
that after a relatively 
short stage of rapid decay driven by the Hall cascade, the 
magnetic field reaches an equilibrium stage which corre¬ 
sponds to the so-called Hall attractor where the Hall cascade 
stops. Later the magnetic field continues to decay mostly 
due to the Ohmic dissipation. So, the field evolution is non¬ 
monotonic with several characteristic time scales in differ¬ 
ent epochs. 

The second possibility is related to existence of two 
regimes of Ohmic decay. In one of them resistivity is due 
to scattering from phonons. This mechanism is important 
till the crust of a NS is hot. In another, resistivity is deter¬ 
mined by impurities in the matter crust. There is a critical 
temperature Tu when magnetic field decay due to phonons 
stops. The time scale of decay on phonons can be ^ 1 Myr 
and even smaller while a NS is hot. For the decay on impu¬ 
rities the time scale depends on the parameter of impurities 
Q, and it can be very long if the parameter is much smaller 
than unity. 

In our approach on every time step we use the formula 
proposed by Aguilera, Pons & Miralles (2008): 


^ 10^^ - 10^® G. Recently it 
Gourgouliatos & Gumming (2013) 


B{t) 


_ Bq X exp(-f/Tohm)_ 

1 + TOhm/THall[l “ exp(-f/rohm)] 


( 1 ) 


In this model two distinct time scales are defined. The first 
one is related to Ohmic decay, Tohm. and the second one — 
to the Hall cascade, rHaii- The equation also can be modified 
to include some minimal value of the field, at which the 
decay is saturated. 
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Let us consider which of two effects play the more 
important role in the crust of normal radio pulsars. The 
timescale of Ohmic decay is; 


■rOhm = 


dTTCrL^ 


( 2 ) 


Where a is local electric conductivity, L is the typical spa¬ 
tial scale of electric currents (can be the local pressure 


height, see Gumming et al. 20041, and c is the speed of light. 
The timescale of the Hall evolution is: 


'THall = 


AirerieL'^ 


(3) 



with rie is local electron density, e is elementary charge, and 
B is local magnetic field. 

As both timescales have the same dependence on the 
spatial scale - r ^ LF', - we can get rid of it if we consider 
the ratio of two time scales: 


T'Ohm 

THall eUeC 


(4) 


If we use the crust properties from 
(|Cumming et al. 2004|l we can simplify Eq. (01): 


T'Ohm 

THall 


= 2.6 X 



(5) 


exp 


4.3 X 10® years 


This equation is valid till the temperature is higher than Ty, 
i.e. while scattering from phonons is important. For fields 
B ~ few xlO^^ G the two time scales are equal for ages 
~few xlO® yrs (see Fig.lH). 


In our study we develop a model which includes all 
these types of decay mainly following prescriptions from 
( Gumming et al. 2004| l. To fix parameters of the crust we 
assume that the Hall time scale is equal to the scale found 
in ( Igoshev & Popov 2014| l for fields ~fewxlO^^ G. To 
take into account thermal evolution of NSs we use cooling 
curves from dShternin et al. 201 1)0 


To understand the role of each types of decay we make 
several runs when one or another of them are switched 
off. In Fig. |2] we present the field evolution only due to 
Ohmic decay via scattering from impurities. This parameter 
is taken to be small: Q = 0.05. This is a realistic assump¬ 
tions for normal field NSs, oppositely for magnetars there is 
a contribution from the part of the crust in the pasta phase, 
where Q is much larger (see IGullon et al. 20141 and refer¬ 
ences therein for details). Then we include scattering from 
phonons, see Fig. [3 For the chosen parameters we see rapid 
decay which stops at ^ 2 x 10® yrs. 

The time scale for the Hall effect depends on the initial 
field value. It can be shorter or longer than the time scale for 
Ohmic decay. In Fig.|4]we compare the time scales for dif¬ 
ferent types of decay. As we see, for realistic parameters the 


Fig. 3 Ohmic magnetic field decay due to scattering from 
phonons and impurities. At the age ~ 2 x 10® years a NS 
cools down to the critical temperature, Ty, when scattering 
from phonons is suppressed. 



True age (years) 

Fig. 4 Typical timescales of magnetic field evolution. 
Thick line on the bottom represents Ohmic decay time scale 
due to phonons. It is truncated when Ty is reached. Dotted 
line on the top corresponds to Ohmic decay due to impuri¬ 
ties. In the middle three curves are shown for the Hall time 
scales for different initial magnetic fields. From top to bot¬ 
tom: Bq = 10^2 G 3 X 1012 G, and lO^® G. 


Hall effect can be as important as Ohmic decay on phonons 
in the time interval of interest ~ (1-4) x 10® yrs. 


For the Hall effect we have also to consider the 
possibility that the Hall attractor develops. Following 
Gourgouliatos & Gumming (2014a) we assume that the 


stage of attractor starts at 3 roHaii- As we do not probe early 
evolution we neglect possible field growth at early stages 
described by Gourgouliatos & Gumming (2015) In Fig. |5] 
we present our main results. 


In Fig. |5]we present field evolution for the case without 
Hall cascade, with Hall cascade but without attractor, and fi¬ 
nally, the one with all effects included. In addition, we plot 
typical ranges of time intervals when a NS age can be es- 


* We than); Peter Shtemin for providing these data and helpful com¬ 
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In our study we assumed that rnaii is compatible with 


the time scale from (Igoshev & Popov 20141. This influ¬ 


enced our estimates of properties of the layer where currents 
are situated. If this basic assumption is wrong, then the re¬ 
sults might be modified significantly. We plan to study it in 
the future. 

Of course, even if we can explore the whole set of pa¬ 
rameters, we are limited by general model assumptions, and 
some (micro)physics is missing. It would be important to 
probe held properties and evolution with observational data 
to put constraints on models and to guide theoretical studies 
in the right direction. 


5 Conclusions 


Fig. 5 Magnetic held decay in three different regimes. 
From top to bottom: no Hall cascade (solid). Hall cascade 
with the Hall attractor (dotted). Hall cascade without attrac¬ 
tor (dot-dashed). For the Hall cascade roHaii = 5 x 10® 
years. The attractor stage is reached at SroHaii- 


timated from the SNR nebula and from kinematic age, and 
the range in which our method can be applied together with 
some other important time intervals and moments. 

We see that for the chosen set of parameters in the range 
considered by Igoshev & Popov (2014) decay might be due 
to joint effect of the Hall cascade and Ohmic decay on 
phonons. Decay is at first slowed down when temperature 
of NS interiors reaches the critical value Ty, and then (out 
of the considered range) when the Hall attractor is reached. 

Note, that if close-by cooling isolated NSs (the 
Magnificent seven) are descendants of magnetars (see, 
Popov et al. 20101, then at ages ~fewx 10® yrs they can be 
at the stage of the Hall attractor. It would be interesting to 
estimate the octupole fields of these objects. 


4 Discussion 


In this paper we presented preliminary results of the analy¬ 
sis of the origin of the magnetic held decay in normal radio 
pulsars found by Igoshev & Popov 2014 We demonstrate 
that the decay with a time scale ^ 4 x 10® yrs and suppres¬ 
sion at ^ 10® yrs (or slightly earlier) can be related to two 
processes. The first one is Ohmic decay due to phonons. 
It has the time scale in the appropriate range, and can be 
stopped when the NS temperature falls down to the critical 
value, T\j. The second one is related to the Hall cascade. 
Fitting the depth at which electric currents are situated we 
can obtain the correct time scale. The held evolution due 
to the Hall effect can reach the stage of attractor at the NS 
age consistent with findings from ( Igoshev & Popov 2014| l. 
However, we think that the second possibility is less proba¬ 
ble. May be the Hall effect also contributes to the field decay 
in normal radio pulsars (see Fig.|5l), but the main contribu¬ 
tion, according to our analysis, is due the Ohmic decay on 
phonons. A more detailed study will be presented elsewhere 


(Igoshev & Popov, in prep. 
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In this section we have to make several comments on 
different assumptions made in our preliminary study. In 
brief, the results depend significantly on parameters as¬ 
sumed. However, the main conclusion is that it is possi¬ 
ble to explain the rate of held decay found in the article by 


Igoshev & Popov 2014 as well as termination of this rapid 


held evolution, in a realistic physical framework. Now let 
us discuss it in more details. 

At first, note that Ty depends on several parameters, and 
so decay on phonons can be terminated earlier or later. Ther¬ 
mal evolution also can vary from the case we use in our 
calculations. This can be either due to differences between 
different NSs (for example, onset of direct URCA process), 
or due to some imperfections of the theoretical model be¬ 
hind the results of calculations we use. As a result, if T\j is 
reached earlier, then in the range of interest decay is mainly 
due to the Hall effect. Oppositely, if T\j is reached later, then 
the Hall effect does not dominate. 
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